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ON  DISCRETE  FAILURE  MODELS 


V.  J.  Padgett  and  John  D.  Spurrier 
University  of  South  Carolina,  Columbia 

Key  Words  -  Discrete  hazard  functions;  Increasing  failure  rate; 
Decreasing  failure  rate;  Discrete  life  data;  Maximus  likelihood  estimation. 


Smeary  4  Conclusions:  *  In  some  situations,  discrete  failure  time 
distributions  are  more  appropriate  to  model  lifetimes  than  continuous  dis- 
trlbutlons.  Very  few  results  for  the  discrete  case  have  been  given  in  the 


literature.  This  paper  provides  three  families  of  discrete  parametric  life* 

(  2»e\ failure  Fate  )  (Prrrfas/nj  FaiJUit-i 

tile  in-fitti^flFR>JlfoFR)J  and  constant 


time  distributions  vhich  are  quite  versati 


failure  rate  models  to  either  uncensored  or  right-censored  life-test  date. 

(  C  hUf) 

The  maximum  likelihood  estimation^of  parameters,  survival  probabilities,  and 
mean  lifetimes  is  investigated,  end  the  KLEs  are  shown  to  be  easily  computed 
by  simple  numerical  methods.  An  example  is  given  for  each  of  the  ss>dels, 
allowing  the  comparison  of  the  proposed  models.  The  example  illustrates  that 
the  discrete  models  presented  provide  a  better  fit  to  discrete  data 

than  previously  proposed  discrete  distributions.  ^ — ' 


1.  INTRODUCTION 


A  great  eaount  of  research  has  been  performed  for  continuous  life 
distributions.  However,  very  little  has  appeared  in  the  literature  for 
discrete  failure  models.  Discrete  failure  data  arise  in  several  common 
situations:  (i)  A  device  can  be  monitored  only  once  per  time  period  (an 
hour,  a  day,  etc.),  and  the  observation  taken  is  the  number  of  time  periods 

euccessfuly  completed  prior  to  the  failure  «f  the  device,  <ii)  A  piece  of _ 

equipment  may  operate  in  cycles,  and  the  experimenter  observes  the  number 
of  cycles  successfully  completed  prior  to  failure.  Many  other  examples  stay 
be  cited. 

In  situations  where  the  observed  data  values  are  very  large,  in  thousands 
of  cycles,  etc.,  a  continuous  distribution  can  be  an  adequate  model  for  the 
discrete  random  variable.  For  example,  the  Birnbaum-Saunders  [3]  distribu¬ 
tion  arises  in  connection  with  fatigue  failure  due  to  bending  of  metal  alloys. 
However,  when  the  observed  values  are  small,  a  continuous  distribution  may  not 
adequately  describe  a  discrete  random  variable,  and  an  appropriate  discrete 
failure  model  is  desirable. 

A  few  results  on  discrete  life  distributions  have  appeared  in  the  litera¬ 
ture.  Kalbfleisch  &  Prentice  [5,  Ch.  2]  briefly  discuss  discrete  failure 
distributions  that  may  be  obtained  from  well-known  continuous  distributions 
such  as  the  exponential  or  Veibull.  Langberg,  Leon,  Lynch  4  Proschan  [5,6] 
examine  properties  of  the  classes  of  discrete  DFR  and  DFRA  distributions. 
Salvia  &  Bollinger  [9]  consider  the  properties  of  discrete  hazard  functions 
and  analogies  with  the  continuous  case.  Barlow,  Bartholomew,  Brenmer  &  Brunk 
(1]  give  a  nonparametric  estimate  of  an  IFR  function  using  isotonic  regression 


techniques*  Their  estimate  is  easily  modified  for  DFR  functions  end  to  allow 
for  censoring. 

The  purpose  of  this  paper  is  to  further  examine  discrete  hatard  functions 
in  a  parametric  setting.  Estimation  for  the  discrete  hazard  functions  given 
by  Salvia  4  Bollinger  [9]  will  be  discussed  for  both  uncensored  and  censored 
lifetime  data  in  Section  2.  In  addition*  some  other  parametric  discrete 
"failure  models  will  be’  proposed  and  estimation  for  them  will  be  studied  in 
Sections  3  and  4.  It  should  be  noted  that  meaningful  discrete  parametric 
models  generally  present  difficult  estimation  problems  and  may  not  be  founded 
on  a  theoretical  basis  [5*  p.  36].  However*  such  models  ere  still  quite 
important  for  discrete  failure  data. 

Notation 


MLE  Maximum  Likelihood  Estimator 

IFR  Increasing  Failure  Bate 

DFR (A)  Decreasing  Failure  Rate  (Average) 

f^  pml{k)  ■  Pr{X  ■  k) 

Sfc  Sf{k)  -  Pr{X  4  k) 


a, a 
F 


E(X) 

L,  En  L 


discrete  hazard*  f^/S^ 

random  sample  of  n  observations  on  a  discrete  random  variable  X 
a  parameter  and  its  MLE 
a  cumulative  distribution  function 
expectation  of  the  random  variable  X 
a  likelihood  function  and  its  natural  logarithm 


2.  THE  MODELS  OF  SALVIA  AND  BOLLINGER 


Salvia  4  Bollinger  [9]  developed  some  results  for  discrete  failure  models 
analogous  to  the  theory  of  the  continuous  case.  Some  of  tbelr  results  follow 


The  discrete  hazard  rate  (or  failure  rate)  is  related  to  the  pmf  by 

\  •  fk/(fk"*-fk+i  +  ...  )  and  fQ  •  h0,  fj  »  hjCl-bg),..., 

fk  -  h^a-hg)  (l-h^)  . . .  (1-h^j) .  Similarly  the  discrete  Sf  is 

$k  »  (l-hjjXl-hj)  ...  (l-h^^j)  for  k»  1,2,...  .  The  mean  life,  if  it  exists, 
m  m>  k-1 

is  E(X)  ■  [  S  •  [  H  (l-h  ).  It  should  be  noted  that  using  Fort 

k-1  *  k-1  j-0  J 

[4,  p.  32-33],  E(X)  exists  if  there  exists  an  e  >  0  and  an  m  such  that 
for  k  >  m,  kfh^/Cl-h^)]  >~1  +  £.  If  for  k  >m,  kti^/d-h^]  S  1,  then 
E(X)  does  not  exist.  Some  other  results  analogous  to  continuous  failure  models 
were  also  given  in  [9]  for  certain  situations.  For  example,  for  a  distribution  with 
IFR,  hjJhjSkjS,,.,  Sk  S  d-hQ)k  and  E(X)  S  (l-h0)/h0.  These  inequalities 
are  reversed  for  DFR  distributions. 

Salvia  &  Bollinger  [9]  also  gave  illustrations  of  discrete  parametric 
failure  models  for  the  constant  failure  rate,  IFR,  and  DFR  cases. 
a.  Constant  Failure  Rate.  For  h^  •  c,  a  constant  0  <  c  S  1  for  all 

If 

k-0,1,2,...,  f^  *  c(l-c)  ,  a  geometric  distribution  which  is  the  discrete 
analog  of  the  exponential  distribution.  Ve  will  give  the  MLE  of  c  for  the 
three  cases  of  a  complete  random  sample  k^,...,kg  from  f^,  a  type  II  censored 
sample, and  a  type  I  censored  sample. 

For  a  complete  random  sample,  the  log-likelihood  is  in  L(c;  k^,...,kQ)  ■ 
n  *  n 

n  in  c  +  l  k.  in(l-c),  which  is  uniquely  maximized  by  c  ■  n/(n  +  J  k.). 

i-1  1  *  *  .  i-1 

Hence,  the  MLE  of  is  ■  (1-0  ,  k*0,l,2 . 

For  a  type  II  censored  sample,  denote  k^  -  minCk^,  ...tk^J.k^  • 

next  smallest  of  (k^,...,kn),  ...,  k^j  “  rth  smallest  of  (k^,...,^), 

where  some  of  the  k^'s  may  be  equal.  If  r+m  items  fail  at  or  before 

k,  m-0,1, . . .,n-r,  the  log-likelihood  is 


t n  L(c;  1^,...,^)  -  (r-Hn)in  c 

+  (  l  4  (n-r-«)Ck<rj +l))ln(l-c) 

which  is  uniquely  maximized  by 

*  r+m 

c  »  (r+m)/[(r+m)  +  J  k^j  +  (n-r-a)(k^rj4l)].  (2.1) 

For  •  type  I  censored  sample*  where  the  life  test_Xs  terminated  after  a _ 

fixed  number  of  ’’cycles"  or  time  units  T,  the  log-likelihood  function  i6 

r 

£n  L(c;  k  ,...,k  )  ■  T  in  f.  +  (n-r)In  S  ,  where  r  is  the  observed  number 
1  n  i-1  k(i)  T  - 

of  failures  by  the  end  of  T  cycles.  The  MLE  c  for  this  case  is  similar  to 

(2.1)  and  is  omitted.  In  fact,  the  expressions  in  in  L  for  arbitrarily 

right-censored  data,  which  includes  type  1  and  type  II  censoring,  have  the 

same  general  form  for  this  or  any  of  the  discrete  models  which  follow. 

b.  An  IFF  Distribution.  For  1^  •  1  -  c/(k+l),  k-0,1,2 . 0  s  c  s  1, 

fk  -  (k-c+l)ck/(k+l)!,  an  IFR  distribution. 

c.  A  DFR  Distribution.  For  h^  ■  c/(k+l),  k*0,l,2,..,,  0  <  c  s  1, 
ffc  -  c(l-c) (2-c) .. .(k-c)/(k+l) I,  a  DFR  distribution. 

Ve  have  found  that  these  IFR  and  DFR  distributions  are  not  sufficiently 
flexible  to  fit  a  vide  variety  of  settings.  For  example,  for  the  IFR  distri- 

A 

bution  the  MLE  c  -  1  if  all  k^s  >  0,  regardless  of  the  values  of  the  k^s. 
For  the  DFR  distribution  E(X)  exists  only  for  c«l.  Generalizations  of  these 
distributions  and  estimation  results  are  given  in  Section  4.1. 

3.  DISCRETE  FAILURE  MODELS  FROM  CONTINUOUS  DISTRIBUTIONS 
lalbfleisch  &  Prentice  [5,  pp.  35-36]  Indicate  that  failure  time  dat*  are 
sometimes  discrete  either  because  of  the  grouping  of  continuous  data  due  to 


imprecise  measurement  or  because  time  ltaalf  is  dlacrctc.  They  alao  glva 

discrete  distributions  baaed  on  grouping  data  from  continuous  distributions 

into  unit  Intervals.  For  example,  if  we  begin  with  the  Veibull  distribution 

weif  (Bt;a)  •  1  -  exp(-(Bt)a),  t  i  0,  then  the  discrete  distribution  obtained 

by  setting  -  weif  (B(k+l);a)  -  weif  (Bk;a)  -  exp(-Bk)°)  -  exp(-($(k+l))a), 

k-0,1,2,...,  is  a  distribution  for  grouped  failure  data.  If  a  *  1,  ffc 

v 

becomes  ffc  -  exp (-Bk) (1-exp (-B))  »  c(l-c)  ,  k-0,1,2,...,  where  c  -  1  -  exp(-8)j 
which  is  the  geometric  model  (a)  of  Salvia  &  Bollinger  [9].  If  a  <  1,  a 
discrete  DFR  distribution  is  obtained,  and  if  a  >  1,  f^  is  an  IFR  distribution. 

More  generally,  it  is  easy  to  show  that  if  a  continuous  IFR(DFR)  distri¬ 
bution  is  used  to  generate  a  discrete  distribution  f^,  then  is  also 

IFR (DFR) .  However,  the  class  of  continuous  distributions  giving  rise  to  discrete 
IFR,  DFR,  or  constant  failure  rate  distributions  is  not  restricted  to  being 
IFR,  DFR,  or  constant  failure  rate,  respectively.  For  example,  the  pdf  g(x)  - 
B  exp(-Bx)  ♦  2  exp(-B(k+Js))  -  2  exp(-B(k+l)),  k  s  x  s  k  +  h,  end  g(x)  -  0, 
k  +  *s  <  x  <  k+1,  k-0,1,2,...  (and  g(x)  -  0,  x  <  0)  does  not  have  constant 
failure  rate,  IFR,  nor  DFR,  yet  ffc  -  exp(-Bk) (1  - exp(-B)) ,  k-0,1,2,,.,,  which 
has  constant  failure  rate,  as  before. 

4 .  SOME  PROPOSED  DISCRETE  PARAMETRIC  FAILURE  MODELS 

Ve  propose  some  more  general  discrete  failure  models  than  those  of  Salvia 
&  Bollinger  19].  FiTBt,  a  generalization  of  their  IFR  and  DFR  models  which  were 
mentioned  in  Section  2  is  given  and  maximum  likelihood  estimation  of  the  para¬ 
meters  is  investigated.  Then  a  discrete  failure  model  which  can  have  increasing, 
decreasing,  or  constant  failure  rate  according  to  the  values  of  one  of  its 
parameters  is  presented,  and  estimation  of  the  parameters  is  considered. 


It  should  be  noted  also  that  some  of  the  well-known  distributions,  in 
addition  to  the  geoaetric  distribution,  way  be  useful  In  modeling  discrete 


life  data.  For  example,  it  say  be  easily  shown  directly,  or  using  the  results 
on  pp.  76-77  of  Barlow  6  Proschan  [2],  that  the  Poisson  distribution  le  an 
IFR  distribution. 

4.1  Generalised  Salvia  end  Bollinger  Models 

_ The  IFR  and  DFR  models  of  Salvia  6  Bollinger  [9]  given  in_Sectlon  2  will 

be  generalized  by  the  addition  of  a  second  parameter  a.  When  a*l  the 
sodels  reduce  to  those  in  Section  2.  When  a-0  the  models  reduce  to  the 
constant  failure  rate  model  of  Section  2. 

The  IFR  Case  (Model  1) 

Let  h  -  1  -  c/(ak+l),  k-0,1,2,...,  (0  s  c  s  1,  0  s  a) .  The  pmf  is 

k  j, 

f  -  (ak+l-e)ck/  II  (ja+l),  k-0,1,2,...,  end  the  Sf  is 
k  J-0 

k  k’1 

S.»  c  /  II  (aj+1),  k-1,2 .  As  O  the  model  approaches  the  Bernoulli 

k  J-o 

model  with  fQ-l-c  end  f^-c.  As  c  -*■  0,  the  model  approaches  the  degenerate 

distribution  with  f Q  - 1 .  The  parameter  c  •  P^fa  >  °)  *nd  o  is  a  shape 

2 

parameter.  The  pmf  is  nonincreasing  in  k  if  c  >  %  and  a  &  (c-1)  /(2c-l) 

or  if  c  S  V  Otherwise,  the  pmf  increases  and  then  decreases  in  k.  Figure  1 

shows  the  pmfs  for  some  values  of  a  end  c. 

Ve  now  consider  the  estimation  of  the  parameters  from  an  observed  random 

sample  k,,...,k  from  f.  .  The  log-likelihood  function  is 
In  k 

n  n  n  k. 

InL  -  I  k.  tnc  +  J  An(ok. +l-c)  -  J  In(jort-l) .  (4.1.1) 

i-1  1  i-1  i-1  J-0 

Three  cases  will  be  considered. 


Cut  1:  a  Known.  If  all  k^  •  0,  then  inL  •  n  in(l-c),  which  gives 


c  •  0.  Next,  if  at  laast  one  k^  >  0,  from  (4.1.1) 

a  n  n 

InL  -  l  k./c  -  l  1/ (ak  + 1  -  c) . 

®c  i-l  x  i-l  1 


(4.1.2) 


Then  — _  inL  <  0  for  all  0  <  c  <  1.  It  follows  from  examining  the  limits 
be 

of  (4.1.2)  as  c  -*■  0  and  c  -*•  1  that  there  is  a  unique  value  c  which 

n  n  A 

maximises  (4.1.1).  If  all  k  >  0  and  ai  7  (l/k.)  /  J  k.»  then  c-1. 

-  1 -  i-l  g  1 - i-l  1  - 

Otherwise,  c  may  be  approximated  by  solving  inL  -  0  for  c  using 

the  Newton-Raphson  procedure. 

Case  ii:  c  Known.  If  all  k  *a  -  0,  then  (4.1.1)  does  not  involve  a. 

-  g 

Bence,  a  is  arbitrary.  If  at  least  one  k^  >  0,  then  inL  -  0  may  or 
may  not  have  a  solution  for  a. 

A  A 

Case  lii:  Both  a  and  c  Unknown.  If  all  k^'s  •  0,  then  c-0  and  a 
is  arbitrary.  Next,  if  nQ  k^e  -  0  and  the  remaining  k^s  »  1,  then 

A*  A 

c  -  (n-nQ)/n  and  o  •  +  •.  Finally,  if  at  least  one  k^  >  1,  it  is  difficult, 
if  not  impossible  to  show  that  unique  MLE's  always  exist.  However,  it  is 
possible  to  find  the  unique  value  of  c  which  maximizes  (4.1.1)  for  fixed  o. 

A  line  eearch  of  the  variable  a  can  then  be  performed  to  obtain  an  approximate 
maximum  for  (4.1.1). 

For  the  cases  of  type  I  and  type  II  censoring,  the  MLE's  for  this  model 
are  obtained  in  the  same  manner.  For  type  I  censoring  after  T  cycles,  the 
log-likelihood  function  is 


r  r 

inL  -  in(c)  l(n-r)  ♦  I  k.,J  +  l  in  [akm  +  1  -  c] 

i-l  W  i-l  w 

r  k,..  T-l 

-  I  lw  in(oj+l)  -  (n-r)  J  in(aj+l). 

i-l  J-0  J-0 


(4.1.3) 


8 


For  type  II  censoring  with  r+m  failures  et  or  before  k^,  ■“0,1 . n-r, 

the  log-likelihood  is 


InL  -  tn(c)[(n-B-r)(k(r)+l)  +  l  k(1))  +  I  tn(ok(1)+l-c) 

(4.1.4) 

k  k 

-  I  f(i)  in(ja+l)  -  (n-m-r)  J(r)  in(joH-l). 
i-1  j-0  3-0 

The  derivatives  of  (4.1.3)  end  (4.1.4)  have  the  same  behavior  as  those  of  (4.1.1), 

so  the  MLE's  may  be  obtained  numerically  as  described  for  uncensored  data. 

The  mean  lifetime  for  this  model  from  the  formula  in  Section  2  is 
•  k-1 

E(X)  -  J  c  /  II  (aj+1)  which  may  be  easily  calculated  to 

k-1  j-0 

any  desired  degree  of  accuracy  by  summing  an  appropriate  number  of  terms  of  the 

A  A 

series.  An  estimate  of  E(X)  is  found  by  replacing  a  end  c  by  a  end  c. 

The  DFR  Case  (Model  2) 

Let  hj^  -  c/(ak+l),  k-0,1,2,...,  (0  <  c  s  1,  0  S  a) .  The  pmf  is  fg-c 


f  -  c  II  (Jor+l-c)/  H  (jaH-1),  k-1, 2,..., 

*  j-0  3-0 

k-1  -1 

and  the  Sf  is  S.  -  II  (oj+l-c)/(aj+l),  where  II  5  1.  As  c  1,  the 

3-0  j-0 

model  approaches  the  degenerate  distribution  with  f  ^  -  1.  The  parameter 

c  -  Pr{X-0}  and  a  is  a  shape  parameter.  The  pmf  is  decreasing  in  k  for 

all  a  and  c.  Figure  2  shows  the  pmf 's  for  various  choices  of  o  and  c. 

Ve  now  consider  the  estimation  of  the  parameters  from  an  observed  random 

sample  k^,...,kg  from  f^.  The  log-likelihood  function  is 

n  k-1  n  k. 

ltd.  -  nine  +  l  l  in(3a+l-c)  -  I  J1  ln(ja+l),  (4.1.5) 

i-1  3-0  i-1  j-0 

-1 

where  J  =  0. 

3-0 


The  results  for  maximising  (4.1.5)  or*  analogous  to  th*  IFR  cat*.  If  o  is 

A  A 

knovn,  there  It  a  unique  0  <  e  £  1  which  maximizes  (4,1.5).  This,  c  ca 

u  a 

he  approx  lasted  by  the  Newton-Raphton  procedure.  If  c  It  known,  then 

a 

^  fnL  •  0  aay  or  may  not  have  a  solution.  If  both  a  end  c  are  unknown, 
then  an  approxinate  maximum  for  (4.1.5)  It  found  by  performing  a  line  search 

#4 

on  a  with  c  -  c  . 

a 

For  the  cases  of  type  I  and  type -II  censoring  the  MLE'-s  are  obtained  in 
the  sane  Banner.  For  type  I  censoring  after  T  cycles,  the  log- livelihood 


function  16 


r  k,.*-l 

T  r  T\*)  , 


£nL  -  rlnc  +  J  I  £n(jc*-l-c)  -  Jwia(JaH)] 

i-1  j-0  j-0 


+  (n-r)  \  JlnKaj+l-c)/(aj+l)l» 

j-o 


(^.1.6) 


For  type  II  censoring  with  r+o  failures  at  or  before  k.  m-0,1 , . ..,n-r. 


the  log-likelihood  is 

InL  -  (srt-r)£nc  +  J  I  J(1)  £n(jO+l-c)  -  J(i)  £n(jo+l)) 
i-1  j-0  j-0 


(4 .1.7) 


+  (n-m-r)  J  r  ia[(aj+l-c)  /  (ttj+1)]. 

J-0 

The  KLEs  may  be  obtained  numerically  as  for  uncensored  data. 

Using  the  results  in  Section  2,  E(X)  exists  if  c-1  or  1  >  c  >  o 
and  does  not  exist  if  c  s  a  and  ci*l.  If  c  >  a,  then  the  mean  lifetime 

may  be  easily  approximated  by  summing  an  appropriate  number  of  terms  of 

•  k-1 

E(X)  -  J  H  (aj+l-c)  /  (aj+1).  An  estimate  of  E(X)  is  found  by  replacing 

k-1  j-0 


4.2  A  Dlicr«t«  "Veibull  Hatard"  Model  (Model  3) 

Let  hy,  ■  1  -  axp(-c(k+l)°),  kB0lll2M..|  (c  >  0,  -  •  <  a  <  •).  Since 

this  Is  similar  to  the  Veibull  Cdf,  we  will  call  this  the  discrete 

Veibull  hazard  model,  or  Model  3.  Thus  the  pmf  is 

f  •  (1-exp (-c (k+1 )°) )  H  exp(-ci°) 

*  i»l 

k 

■  exp(-c  J  i°)  (1-exp (-c(k+l)  )),  k»0,l,2,..., 

_  -  i-1  -  - - 

k  0 

and  the  Sf  is  S.  •  exp(-c  2  *)»  k-1,2,...,  where  2=0* 

*  i-1  i-1 

Mote  that  for  a-0,  this  model  reduces  to  the  constant  failure  rate  model 

of  Section  2.  For  a  >  0,  the  distribution  has  IFR,  and  for  a  <  0,  it  has 

DFR.  Thus,  this  discrete  nodel  is  quite  flexible  with  respect  to  choice  of 

failure  rate,  analogous  to  the  Veibull  distribution  in  the  continuous  case. 

The  parameter  a  is  essentially  a  shape  paraaeter  and  c  is  a  "spread" 

parameter.  For  small  values  of  c,  the  probability  is  spread  out  in  the  tail 

of  the  distribution,  whereas  for  large  c,  most  of  the  probability  is  at  k-0. 

For  example,  for  c»5,  fQ  “  0.9933.  For  a  <  0,  the  pmf  is  always  decreasing 

in  k.  For  a  >  0,  the  pmf  is  either  decreasing  or  first  Increasing  and  then 

decreasing  in  k.  In  particular,  for  small  c,  the  largest  ffc  is  at  large 

values  of  k.  For  example,  for  a  -  0.5,  c  *  .01,  the  largest  value  of  f^ 

is  at  k-13.  Figure  3  shows  the  pmf 's  for  some  values  of  a  and  c.  It  appears 

that  for  various  values  of  c  >  0,  a  in  the  range  -1  S  o  s  1  is  sufficient 

to  describe  many  IFR,  DFR,  or  constant  failure  rate  discrete  distributions. 

Ve  now  consider  the  estimation  of  the  parameters  from  an  observed  random 

sample  k^, . . . ,kQ  from  f^.  The  log-likelihood  function  is 

n  n  k. 

AnL  -  2  An[l-exp(-c(k.+l)a]  -  c  2  I  J  • 

*  4-1  4-1 


(4.2.1) 


Three  cese*  will  be  considered 


Case  i:  o  Known.  If  et  least  one  k^  >  0,  from  (4.2.1) 


-r?  InL  -  l  (k.+l)°/  [exp(c(k  +1)°  -  1)  -  l  I1  j°. 
Sc  i-1  x  1  1-1  j-1 


(4.2.2) 


Then,  Ini  <  0  for  all  c  >  0,  and  as  e  0,  (4.2.2)  approaches  +  “, 

3c 

whereas  as  c  •*•  *,  (4.2.2)  approaches  a  negative  constant.  Hence,  there  Is 

g 

a  unique  MLE  for  c.  This  KLE  can  be  found  by  aolving  InL  •  0  numerically 
using  Nevton-Raphson  Iteration. 


If  k^  -  0  for  all  i-l,...,n,  then  £n  L  does  not  involve  a  and  is 


an  Increasing  function  of  c.  Hence,  c  -  +  •  and  a  is  arbitrary.  These 

A  A 

estimates  correspond  to  f^  -  1  and  f^  -  0,  k  >  0. 

Case  ii:  c  Known.  For  any  fixed  value  of  c,  it  can  be  shovn  that  there 

3 

exists  at  least  one  solution  to  InL  -  0  for  a  if  at  least  one  k^  >  1. 

The  solution  may  not  be  unique.  If  all  k^  -  0,  InL  does  not  Involve  a, 
so  a  is  arbitrary.  Finally,  if  all  k^*s  •  1  or  some  k^'s  »  0  and  the 
remaining  k^'s  -  1,  then  there  nay  be  no  finite  value  for  a. 

Case  Hi:  Both  a  and  c  Unknown.  If  at  least  one  k^  >  0,  it  is  difficult 
if  not  impossible,  to  show  that  unique  MLEs  always  exist.  However,  it  is 
possible  to  find  the  unique  values  of  c  which  maximize  (4.2.1)  for  fixed 
values  of  a  using  the  Nevton-Raphson  procedure  for  solving  InL  -  0. 

A  line  search  on  the  variable  o  can  then  be  performed  to  obtain  an  approximate 
maximum  for  (4.2.1). 

For  the  cases  of  type  I  and  type  II  censoring,  the  MLEs  for  this  failure 
model  are  obtained  in  the  same  manner  as  just  described.  For  type  I  censoring 
after  T  cycles,  the  log-likelihood  function  is 


r  k 


ioL  -  l  *nU-«p(-c)<k, .,♦!)“)]  -  c  l  l*'  jl 


i-1  J-l 


-  c(n-r)  J  J°. 

J-l 


(4.2.3) 


For  type  II  censoring  with  r-ta  failures  et  or  before  k^j,  *-0,1,. ..,n-r, 
the  log-likelihood  is 


r+m  „  r-hi  k. .. 

inL  -  l  in[l-axp(-c(k,  .+1)°)]  -  c  J  _JUJ  ja 


i-1  -j-l 


k,  *+l  ^ 
-  c(n-r-m)  ja. 

J-l 


(4.2.4) 


The  derivatives  of  (4.2.3)  end  (4.2.4)  have  the  same  behavior  a6  those  of 

(4. 2.1), so  the  MLEs  may  be  obtained  numerically  as  described  for  uncensored  data 

•  k 

The  mean  lifetime  for  this  model  is  E(X)  -  J  exp(-c  J  jj .  By  the 

k-1  J-l 

results  referenced  in  Section  2 ,  E(X)  exists  for  ell  c  >  0  end  a  >  -1 
end  does  not  exist  when  a  <  -1  for  any  value  of  c  >  0.  For  a  -  -1,  E(X) 
exists  when  c  >  1  but  does  not  exist  for  c  £  1.  When  E(X)  exists,  it  may 
be  easily  obtained  to  any  desired  degree  of  accuracy  by  summing  an  appropriate 
number  of  terms  in  the  series.  An  estimate  of  E(X)  can  be  found  by  using 

A  * 

o  end  c.  For  o  <  0  (DFR),  a  lower  bound  for  E(X)  is  exp (-c)/ [1-exp (-c)J « 

* 

which  can  be  estimated  by  replacing  c  with  c.  For  a  >  0  (IFR),  this  quantity 


is  an  upper  bound  for  E(X). 


5.  EXAMPLES 


Kelson  [8]  reported  times  to  breakdown  in  minutes  of  an  insulating  fluid 
subjected  to  various  test  voltages.  We  took  as  an  example  the  failure  times 
at  36  hr  end  reduced  each  observation  to  the  next  lowest  integer,  which  corres¬ 
ponds  to  observing  the  life  test  once  per  minute  and  recording  the  lifetime  as 
*V  if  the  breakdown  occurred  in  time  interval  fk,k+l).  Kelson  (8]  assumed 


an  underlying  Veibull  distribution  for  the  continuous  fellure  Claes.  As 
examples,  we  will  consider  estiastlon  from  this  dete  for  the  "discrete  Veibull" 
aodel  of  Section  3,  Models  2  end  3  of  Section  4,  end  the  nonpereaetric  DFR 
aodel  of  Barlow,  Bartholomew,  Br earner,  end  Brunk  [1]. 

The  dete  that  we  use,  es  edepted  from  Nelson  [8]  for  36kV,  ere: 
0,0,0,0,1,1,2,2,2,2,3,3,5,13,25  (minutes  before  fellure).  The  coaputed  MLEe  of 
peremeters,  pmfs,  Sfs,  discrete  hezerds,  end  ween  lifetimes  for  the  various 
models  ere  shown  for  comparison  In  Table  1.  A  hazard  plot  of  these  15  observa¬ 
tions  indicated  that  the  life  distribution  is  DFR.  This  is  supported  by  the 

* 

MLE  of  a  in  the  "discrete  Veibull  hazard"  model  (Model  3)  which  was  o  c  -.39: 

It  can  be  seen  from  Table  1  that  the  nonparametric  DFR  model  yields  the 
largest  value  for  £nL.  This  is  to  be  expected  since  maximization  of  £nL 
is  done  over  all  DFR  distributions  (Including  the  other  models  in  Table  1). 

The  nonparametric  model  has  two  drawbacks,  particularly  when  the  number  of 

A 

failures  is  small.  First,  the  estimate  of  the  hazard  function,  h^,  decreases 
at  k  only  if  there  is  an  observed  failure  at  k.  Secondly,  for  k  larger  . 
than  the  largest  failure  tine,  h^  Is  not  uniquely  defined.  For  the  purposes 
of  Table  1  ^  =  h2J  for  k  >  25. 

The  three  parametric  models  yield  similar  estimates  of  the  failure  distri¬ 
bution.  Ve  make  the  following  distinctions  among  the  models.  In  terms  of 
£nL,  Model  2  performs  the  best  for  this  data.  Calculation  of  the  estimators 
is  more  complicated  for  the  discrete  Veibull  model  than  for  Models  2  and  3. 

Model  3  and  the  discrete  Veibull  aodel  are  more  flexible  than  Model  2  since  they 
allow  both  IFR  and  DFR  distributions.  These  proposed  models  provide  a  large 
class  of  distributions  for  fitting  discrete  failure  data. 
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